Strong continuous- variable entanglement in modulated quantum dissipative dynamics 
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We investigate the creation of entangled states of bright light beams obeying the condition of 
strong Einstein-Podolsky-Rosen-like paradox criterion in time-modulated quantum dissipative dy- 
namics. Having in view the generation of these states we propose a non-degenerate optical para- 
metric oscillator(NOPO) driven by an amplitude-modulated pump field. We develop semi-classical 
and quantum theories of this device for all operational regimes concluding that, contrary to ordi- 
nary NOPO, the continuous- variable entanglement becomes approximately perfect for high level of 
modulation. Our analytical results are in well agreement with numerical simulations. 
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Continuous-variable (CV) entangled states of light 
beams provide excellent tools for testing the foundations 
of quantum physics and arouse growing interest due to 
apparent usefulness as a promising technology in quan- 
tum information and communication protocols 0, 0- 
The efficiency of quantum information schemes and quan- 
tum measurements significantly depends on the degree of 
entanglement of used states. On the other hand, in the 
majority of real applications bright stable light beams 
are required. It is therefore highly desirable to elabo- 
rate reliable sources of light beams having the mentioned 
properties. 

The recent development of CV quantum information is 
stipulated mainly by preparation of two-mode squeezed 
states, which are an approximation to EPR (Einstein- 
Podolsky-Rosen) entangled states and can be easilygen- 
erated by a nondegenerate parametric amplifier [J 0- 
However, up to now the generation of bright light beams 
with high degree of CV entanglement meets serious prob- 
lems. One of these is the degradation of entanglement 
due to uncontrolled dissipation and decoherence pro- 
cesses. Nevertheless, an experimental progress in the 
generation of CV entangled intensive light beams from 
NOPO has been reported in 5]. A source of CV en- 
tangled light has been also built using two quadrature 
squeezed beams combined on a beam splitter |2( , and ex- 
perimental characterization of CV entanglement has been 
demonstrated using this method in |(j . CV entanglement 
of phase-locked light beams was recently studied in Q • 

Note that because the number of experimentally ac- 
cessible and feasible multi-wave nonlinear interactions 
leading to formation of entanglement is rather limited, 
the class of schemes that may be practically elaborated 
is restricted. In this Letter we point out that the 
class of currently proposed schemes for generation of in- 
tensive light bemas with high degree of entanglement 
may be significantly extended if instead of monochro- 
matic pumping wc consider periodically modulated pump 
fields. As a realization of this idea, we propose in this 
Letter a novel scheme of NOPO in a cavity, driven by 



an amplitude-modulated electromagnetic field and stress 
that this scheme provides highly effective mechanism for 
improvement of the degree of CV entanglement, even in 
the presence of dissipation and cavity induced feedback. 

CV entangling resources are usually analyzed based on 
two-mode squeezing through the variances of the quadra- 
ture amplitudes of two generated modes. In NOPO, the 
two-mode integral squeezing, which characterizes the CV 
entanglement, reaches only 50% relative to the level of 
vacuum fluctuations, if the pump field intensity is close 
to the generation threshold ja, li| • We will demonstrate 
here, that the level of two-mode squeezing in the pro- 
posed scheme is not limited, which indicates a high de- 
gree of quadrature entanglement obeying the condition 
of EPR-like paradox criterion as quantified by Reid and 
Drummond [2j. 

It seems intuitively clear that such an achievement 
is due to the control of quantum dissipative dynam- 
ics through the application of suitably tailored, time- 
modulated driving field. Some examples of the suppres- 
sion of quantum decoherence as well as an improvement 
of nonclassical statistics of oscillatory excitation numbers 
have been considered in . 

We consider a type-II phase-matched NOPO with 
triply resonant cavity that supports the pump mode and 
two orthogonally polarized modes of subharmonics. The 
Hamiltonian describing intracavity interactions within 
the framework of rotating wave approximation and in 
the interaction picture is 

H = ihf(t)(e i( -^-^at-e- i( -^- u ^a^ 

+ihk (e %i5>h a^a^ a% - e^ k a^aia 2 ) , (1) 

where a, are the boson operators for cavity modes at 
the frequencies Wj. The pump mode 03 is driven by 
an amplitude-modulated external field at the frequency 
u)l = W3 with amplitude f(t) that is a periodic function 
with modulation frequency iCwj,. The constant fce l * fc 
determines an efficiency of the down-conversion process 
in %( 2 ) medium. 



The system of interest is dissipative because the sub- 
harmonic modes suffer losses. Taking into account the 
cavity damping rates 7^ of the modes we consider the 
case of zero detunings and high cavity losses for the pump 
mode (73 3> 7) under the assumption that 71 = 72 = 7- 
However, in our analysis we allow for the pump deple- 
tion effects. Following the standard procedure of quan- 
tum optics we derive in the positive P-repesentation [Jjj 
the stochastic equations for two groups of independent 
complex c-number variables cti. 2 and /?i >2 corresponding 
to operators ai j2 and a^ 2 



^ = -( 7 + \a 2 p 2 ) ai + e (t) f3 2 + W ai (t) , (2) 
^ = -(7 + Aa 2/ 3 2 )/?i + e* (t) a 2 + W 01 (t) . (3) 



Here: e(t) — f(t)k/-f 3 , A = fc 2 /73 aim equations for 
«2, 02 are obtained from (J2J, © by exchanging the sub- 
scripts (1)^(2). Our derivation is based on the Ito 
stochastic calculus, and the nonzero stochastic correla- 
tions are: (W ai (t) W a2 (*')) = (e (*) - Ac*ia 2 ) S (t - t'), 
{W 01 (t) W(} 2 (if)) = (e (i) - \/3if3 2 ) 8(t- t'). Note, that 
while obtaining these equations we used the transformed 
boson operators 04 — > aiexp (— with $i being $3 = 
$l, *i = $2 = I ($l + $fc)- This leads to cancella- 
tion of phases at intermediate stages of calculation. As 
a result, the equations depend only on real and positive 
coupling constants. 

The equations of motion are with time-dependent co- 
efficient. Nevertheless, surprisingly, it is possible to find 
their analytical solution in the semiclassical approach for 
an arbitrary, but real modulation amplitude f(t). The 
presentation of both the semiclassical and quantum the- 
ories of such a time-modulated NOPO is another impor- 
tant goal of this Letter. 

First, we shall study the solution of stochastic equa- 
tions in semiclassical treatment, neglecting the noise 
terms, for mean photon numbers rij and phases tpj of the 
modes (rij — aj(3j, ipj — ln(ctj//3j)) for time-intervals 
exceeding the transient time, t 3> 7 _1 - An analysis 
shows that similar to the standard NOPO, the con- 
sidered system also exhibits threshold behavior, which 
is easily described through the period-averaged pump 
field amplitude J(tj = ± £ f(t)dt, where T = 2n/S. 
The below-threshold regime with a stable trivial zero- 
amplitude solution is realized for / < f t h, where fth = 
773/fc is the threshold value. When / > fth, the sta- 
ble nontrivial solution exists with the following proper- 
ties. First, the mean photon numbers for subharmonic 
modes n i = (a^ai) — \ai\ 2 are equal one to the other 
(noi = ?io2 = no) due to the symmetry of the system, 
7i = 72 = 7- As for usual NOPO, the phase difference 
is undefined due to the phase diffusion, while the sum of 
phases is equal <*p\-\-<*pi = 2^k. The mean photon number 
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FIG. 1: The semiclassical mean photon number versus dimen- 
sionless time for the parameters: fc/7 — 5 ■ 1CP 4 , 73/7 = 25, 
S/j = 2, J = 3/ th : fx = (curve 1), /1 = 0.4/ (curve 2) and 
fx = 1.2/ (curve 3). 



no(t) satisfies the following equation 



—n (t) = 2n (t) (e (t) 



7 - An (t)) , 



(4) 



which can be explicitly integrated, giving a periodic so- 
lution 

n \t) = 2\J exp(^2 J (e (t 1 + 1) - 7) dt^j dr. (5) 

Note, that in J5| the terms corresponding to the tran- 
sient dynamics are omitted for simplicity. In the ab- 
sence of modulation rio(t) reaches the steady-state n — 
(/ — fth)/k in accordance with the theory of standard 
NOPO for zero detunigs. 

We emphasize that this model is available for ex- 
periments. It can be implemented at least for NOPO 
driven by polychromatic pump fields. Particularly, for 
the following pump field E ext {t) = f cos(wz,i + $l) + 
(/i/2) (cos((w L + S)t + $ L + $) + cos((w L - S)t + 
the system is described by the Hamiltonian with har- 
monic modulation amplitude f(t) = f + fiCOs(St + $). 
We present below the final results for the case of such 
harmonic modulation assuming without loss of general- 
ity / > 0, /1 > and $ = 0. In this case the mean 
photon number JSJ reads as 



^ 1 (i) = 2 A f cxp ^2 7 t 



J_ 

fth 



exp 



[sin (5 (t + r)) - sin {St)]) dr. (6) 

V ° fth J 



This result is illustrated in Fig.Q for the different levels 
of modulation. 

To characterize the CV entanglement we address to 
both the inseparability criterion proposed by Duan et 
al. and Simon and the EPR paradox criterion 0|. 
These criteria could be quantified by analyzing the vari- 
ances of a relevant distance V- — V (X\ — X 2 ) and to- 
tal momentum V+ = V(Ji+i^) in the terms of the 



3 



quadrature amplitudes of two modes X k — X k (Ofc) = 
-^(a+e-* 6 * + a k e ie "), Y k = X k (Q k - § ) , (k = 1,2), 

where V(x) — (a; 2 ) — (x) 2 is a denotation of the vari- 
ance. The inseparability criterion, or weak entanglement 
criterion, for the sum of variances quite generally reads 
as V+ + V- < 2, and due to the mentioned symmetries is 
reduced to the following form V = V+ = V- < 1, while 
for the product of variances this criterion has the form 
V+V- = V 2 < 1. The strong CV entanglement criterion 
shows that when the inequality V+V- < 1/4 is satisfied, 
there arises an EPR-like paradox. Obviously, the suf- 
ficiency condition for inseparability is weaker than the 
EPR condition. Both criteria have been used to charac- 
terize the entanglement mainly in spectral measurements 
0. Contrary to that, we confine ourselves to analyzing 
only the total intracavity variances. These variances are 
expressed through the stochastic variables using the re- 
lationships between normally-ordered operator averages 
and stochastic moments with respect to the P-function. 
Restoring the previous phase structure of intracavity in- 
teraction, we obtain that V+ = V- = V and 



V=l+ (aift) + (a 2 (3 2 ) - (ai« 2 ) e 4W - e~ lW ,(7) 

where 6 = 6i + 6 2 + <&£, + $fc- 

Further, we will calculate the variance in the standard 
linear treatment of quantum fluctuations. To this end, it 
is convenient to use the following moments of quantum 
variables, which are given below for both operator- and 
P-representations: (ri + ) := (a^ ai) + (a^a^) = (otxPx) + 
(a 2 f3 2 ), (R) ■= (at a l) + ( a t a 2) - (a\a 2 ) - (a+a^) = 
((ai-ftJGSi-aa)), (Z) := ((a+ ai - a+a 2 ) 2 ) = 

(^(aiPi - a 2 (3 2 ) 2 ^ + (aiPi) + (a 2 f3 2 ). As can be seen, 
the possible minimal level of variance, realized under ap- 
propriate selection of phases 81 + 62 = — $l — ®k in 
formula J7J, is expressed through the moment (R) as 
V(t) = 1 + (R(t)). Using Ito rules for changing the 
stochastic variables, we obtain from © , © the following 
equations: 

|(n+) = {2e(t)-2 1 -X)(n+)-X(n 2 + ) 

- 2s (t) (R) + X (Z) , (8) 
= - (2e (t) + 2 7 + A) (R) - A (n+R) 

- 2e(t) + X(Z), (9) 
= -4 7 (Z) + 2 7 (n+). (10) 



> 



From Ea. lfTrHl (Z) can be expressed as a function of (n+). 
Substituting this expression into ©, © we get the fol- 
lowing equations which are convenient for the perturba- 
tive analysis of quantum fluctuations 



dt 



(n+) = (2s{t)-2 1 -X)(n + )-X(nl)-2s(t)(R) 



> 



2 7 A / e^ T -^(n + (T))dT, 

J — OO 

- (2e (t) + 2 7 + A) (R) - A (n+R) 
2e (t) + 2 7 A / e^ T ~V (n+ (r)) dr. 



(11) 



(12) 



First, we consider the above-threshold regime lineariz- 
ing quantum fluctuations around the stable scmiclassical 
solutions. In the linear treatment of quantum fluctua- 
tions we have the expansions (n+) = nio + n 2 Q + (Sn+) = 
2n + (Sn+), (R) = R° + (SR) = (SR), {n+R) = 2n Q (SR), 
(n + ) = 4no (Sn+), where it was assumed that nio = 
n 2 o = no(t), if 1 + ip 2 = 2nk, and hence R° = 0. Note, 
that in the current experiments the ratio of nonlinear- 
ity to dumping is small, fc/7 1 (typically 10~ 6 or less), 
and hence A/7 = k 2 / (773) -C 1 is the small parameter of 
the theory. Therefore, the zero order terms in the above 
expansion correspond to a large classical field of the or- 
der 7/ A in accordance with Eq.©, while the next terms 
describing the quantum fluctuations are of the order of 1. 
On the whole, combining the procedure of linearization 
with A/7 -C 1 approximation we get a linear equation for 
the variance V(t) = 1 + (SR) 



dt 



V(t) = -2( 1 + e(t) + Xn Q (t))V(t) + 2Xn Q (t) 

+ 2 7 + 4 7 A / e^-^noMdr, (13) 



with the following periodic asymptotic solution 

V (t) — 2 J exp (-2 J (7 + e (*') + Xn (t')) dA x 



7 + Xn (r) + 27A 



°M T '- T )n Q (T')dT' 



dilA) 



It should be noted, that the result l|14|) is also obtained 
from initial equations ©, © in the photon number and 
phase variables using, however, more complicated calcu- 
lations. In the absence of modulation the formula (|14|) 
coincides with an analogous one for ordinary NOPO. The 
analysis of the below-threshold regime is more simple and 
leads to formula (|14f> with no = 0. 

Let us outline some conclusions from this result for 
the case of harmonic modulation, using that e (t) = 
^ {f + f 1 cos (St)) and the expression ©. Typical re- 
sults are presented in Fig.© for the above-threshold 
regime, f / fth — 3. The variance is seen to show a time- 
dependent modulation with a period 2-k/8. The dras- 
tic difference between the degree of two-mode squeez- 
ing/entanglement for modulated and stationary dynam- 
ics is also clearly seen in Fig.©. The stationary variance 
(curve 1) near the threshold having a limiting squeezing 
of 0.5 (see also Fig.©, curve 1) is bounded by quan- 
tum inseparability criterion V < 1, while the variance 
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FIG. 2: The variance V(£) given by the linear theory, versus 
the dimensionless time ■yt for the parameters as in Fig.Q. 
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FIG. 3: The minimum level of the variance versus f /fth for 
three levels of modulation: fi ~ (curve 1), /i = 0.75/ 
(curve 2) and /i = 2/ (curve 3). The parameters are: k/f = 
5 ■ 1CT 4 , 73 / 7 = 25, 8h = 2. 



for the case of modulated dynamics obeys the EPR cri- 
terion V 2 < 1/4 of strong CV entanglement for definite 
time intervals. In particular, the minimum values of the 
variance and corresponding photon numbers of Fig. Q at 
fixed time intervals t m = to + 2Trm/S, (m = 0, 1, 2...) are: 
n ~ 6.16 • 10 7 , V mm ~ 0.27, t = 2.647" 1 (curve 3) and 
n ~ 1.71 ■ 10 8 , V mm ~ 0.56, t = 2.51-f" 1 (curve 2). 
The dependence of V m in on the period-averaged pump 
field amplitude is shown in Fig.© for different levels of 
modulation. As it is expected, the degree of EPR entan- 
glement increases with ratio fi/f. Another peculiarity 
here is that the stationary variance (curve 1) has a char- 
acteristic threshold behavior (see, for example 8]) that 
disappears in case of strong modulation (curve 3). Wc 
conclude, that here nothing of the kind of principle is 
valid that prohibits the reaching of approximately per- 
fect squeezing/entanglement for the case of strong mod- 
ulation. The production of strong entanglement occurs 
for the period of modulation comparable with the char- 
acteristic time of dissipation, 5 ~ 7. For both asymptotic 
cases of slow (S <C 7) and fast (S ^> 7) modulations this 
effect disappears. 

The linearized theory is applicable only out- 
side the critical region. As our analysis shows, 
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FIG. 4: Comparison of the squeezed variance between the 
analytical solution (curve 1) and numerical simulation (curve 
2). The parameters are: / = f t h, fi = 0.5/, A/7 = 0.01, 
73 /7 = 25,<5/7 = 2. 

the condition of the validity of linear results for 
the near-threshold regimes reads as \f / fth — 1 
(A/7)exp[2(/i//th)(7/<y)]. For typical A/7 « 1, this 
condition is fairly easy to satisfy even for narrow criti- 
cal ranges provided that 5 » 7. As a rule, the quantum 
corrections diverge at the classical threshold, although, 
the variance l|14fl surprisingly is well defined also at the 
threshold. Nevertheless, in order to verify the accuracy of 
our analytical calculations we performed numerical simu- 
lation based on the quantum state diffusion method (|l3|. 
for the case of time-modulated fields, see ^3)- This sim- 
ulation tends to disagree with our analytical results for 
comparatively high values of A/7 parameter, especially at 
the threshold, as is shown in Fig. 10}. Generally speak- 
ing the divergence between the numerical and analytical 
results is typically in the order of A/7. 

In conclusion, we have proposed a new approach to 
generation of strongly EPR entangled states of intense 
light beams based on the time-modulation of quantum 
dissipative dynamics. We have demonstrated an essen- 
tial improvement of the degree of entanglement in mod- 
ulated NOPO in comparison with the ordinary one, if 
the frequency of modulation is close to the decay rate 
of dissipative precesses. We have shown that surpris- 
ingly the CV entanglement becomes approximately per- 
fect for high-level of modulation, but is only realized for 
a periodic sequence of time intervals synchronized with 
the modulation period. The analytically predicted re- 
sults were confirmed in our numerical simulation. In our 
study we have not analyzed in detail any physical mech- 
anisms leading to strong EPR entanglement due to the 
modulation. This topic is currently being explored and 
will be the subject of forthcoming work. We believe that 
the results obtained are applicable to a general class of 
quantum dissipative systems and can serve as a guide for 
further theoretical and experimental studies of intensive 
light beams with high degree of entanglement. 
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